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The goal of ultra-rclativistic nuclear collisions is the 
creation and study of the quark-gluon plasma (QGP), 
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a state of matter whose existence at high energy den- 
sity is predicted by Quantum Chromodynamics. One 
of the experimental observables that is sensitive to the 
properties of this matter is the azimuthal distribution 
of particles in the plane perpendicular to the beam di- 
rection. When nuclei collide at finite impact parameter 
(non-central collisions), the geometrical overlap region 
and therefore the initial matter distribution is anisotropic 
(almond shaped). If the matter is interacting, this spa- 
tial asymmetry is converted via multiple collisions into an 
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anisotropic momentum distribution [T]. The second mo- 
ment of the final state hadron azimuthal distribution is 
called elliptic flow; it is a response of the dense system to 
the initial conditions and therefore sensitive to the early 
and hot, strongly interacting phase of the evolution. 

At RHIC large elliptic flow has been observed and is 
one of the key experimental discoveries [2HS] ■ Theoretical 
models, based on ideal relativistic hydrodynamics with a 
QGP equation of state and zero shear viscosity, fail to 
describe elliptic flow measurements at lower energies but 
describe RHIC data reasonably well [TJ. Theoretical ar- 
guments, based on the AdS/CFT conjecture [5], suggest 
a universal lower bound of 1/Att |J| for the ratio of shear 
viscosity to entropy density. Recent model studies incor- 
porating viscous corrections indicate that the shear vis- 
cosity at RHIC is within a factor of ~5 of this bound [TuT - 

The pure hydrodynamic models [7J 1141 115) and mod- 
els which combine hydrodynamics with a hadron cascade 
afterburner (hybrid models) [TBI HI] that successfully de- 
scribe flow at RHIC predict an increase of the elliptic flow 
at the LHC ranging from 10% to 30%, with the largest 
increase predicted by models which account for viscous 
corrections jTSTUB] at RHIC energies. In models with vis- 
cous corrections, t>2 at RHIC is below the ideal hydrody- 
namic limit |12l 117] and therefore can show a stronger in- 
crease with energy. In hydrodynamic models the charged 
particle elliptic flow as a function of transverse momen- 
tum does not change significantly [7J H3] , while the pt- 
integratcd elliptic flow increases due to the rise in average 
Pt expected from larger radial (azimuthally symmetric) 
flow. The larger radial flow also leads to a decrease of 
the elliptic flow at low transverse momentum, which is 
most pronounced for heavy particles. Models based on a 
parton cascade [19 , including models that take into ac- 
count quark recombination for particle production [20j . 
predict a stronger decrease of the elliptic flow as func- 
tion of transverse momentum compared to RHIC ener- 
gies. Phenomenological extrapolations [3T] and models 
based on final state interactions [53] that have been tuned 
to describe the RHIC data, predict an increase of the el- 
liptic flow of ~ 50%, larger than other models. A mea- 
surement of elliptic flow at the LHC is therefore crucial 
to test the validity of a hydrodynamic description of the 
medium and to measure its thermodynamic properties, in 
particular shear viscosity and the equation of state [23] . 

The azimuthal dependence of the particle yield can be 
written in the form of a Fourier series [23J [35] : 

d 3 N 1 d 2 N ( ^ \ 

where E is the energy of the particle, p the momentum, p t 
the transverse momentum, <f> the azimuthal angle, y the 
rapidity, and the reaction plane angle. The reaction 
plane is the plane defined by the beam axis z and the 
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impact parameter direction. In general the coefficients 
v n = (cos[n((f> — ~*&r)]) are pt and y dependent - therefore 
we refer to them as differential flow. The integrated flow 
is defined as an average evaluated with d 2 N/dp t dy used 
as a weight. The first coefficient, vi, is called directed 
flow, and second coefficient, V2, is called elliptic flow. 

We report the first measurement of elliptic flow of 
charged particles in Pb-Pb collisions at the center of 
mass energy per nucleon pair ^/s NN = 2.76 TeV, with 
the ALICE detector [25H28] . The data were recorded in 
November 2010 during the first run with heavy ions at 
the LHC. 

For this analysis the ALICE Inner Tracking System 
(ITS) and the Time Projection Chamber (TPC) were 
used to reconstruct the charged particle tracks. The 
VZERO counters and the Silicon Pixel Detector (SPD) 
were used for the trigger. The VZERO counters are 
two scintillator arrays providing both amplitude and 
timing information, covering the pseudorapidity range 
2.8 < n < 5.1 (VZERO-A) and -3.7 < n < -1.7 
(VZERO-C). The VZERO time resolution is better than 
1 ns. The SPD is the innermost part of the ITS, con- 
sisting of two cylindrical layers of hybrid silicon pixel as- 
semblies covering the range of \i]\ < 2.0 and \r]\ < 1.4 for 
the inner and outer layer, respectively. The minimum- 
bias interaction trigger required at least two out of the 
following three conditions [25]: i) two pixel chips hit in 
the outer layer of the silicon pixel detectors ii) a signal in 
VZERO-A iii) a signal in VZERO-C. The bunch inten- 
sity was typically 10 7 Pb ions per bunch and each beam 
had 4 colliding bunches. The estimated luminosity was 
5 x 10 23 cm -2 s _1 , producing collisions with a minimum 
bias trigger at a rate of 50 Hz including about 4 Hz nu- 
clear interactions, 45 Hz electromagnetic processes and 1 
Hz beam background. 

A removal of background events was carried out of- 
fline using the VZERO timing information and the re- 
quirement of two tracks in the central detector. A study 
based on Glauber model fits to the multiplicity distribu- 
tion (see also [25J) in the region corresponding to 80% of 
the most central collisions, where the vertex reconstruc- 
tion is fully efficient, allows for the determination of the 
cross section percentile. Only events with a vertex found 
in \z\ < 10 cm were used in this analysis to ensure a 
uniform acceptance in the central pseudorapidity region 
\n\ < 0.8. An event sample of 45 x 10 3 Pb-Pb collisions 
passed the selection criteria and was used in this analysis. 
The data are analyzed in centrality classes determined by 
cuts on the uncorrected charged particle multiplicity, in 
pseudorapidity acceptance \n\ < 0.8. Figure [I] shows the 
uncorrected charged particle multiplicity in the TPC for 
these events and indicates the nine centrality bins used 
in the analysis. 

To select charged particles with high efficiency and to 
minimize the contribution from photon conversions and 
secondary charged particles produced in the detector ma- 
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FIG. 1. The uncorrected multiplicity distribution of charged 
particles in the TPC (\w\ < 0.8). The centrality bins used 
in the analysis are shown and the cumulative fraction of the 
total events is indicated in percent. The bins 60-70% and 
70-80% are not labeled. 



terial, the following track requirements were applied for 
tracks measured with the ITS and TPC. The tracks are 
required to have at least 70 reconstructed space points 
out of the maximum 159 in the TPC and a (x 2 ) per 
TPC cluster < 4 (with two degrees of freedom per clus- 
ter) . Additionally, at least two of the six ITS layers must 
have a hit associated with the track. Tracks are rejected 
if their distance of closest approach to the primary vertex 
is larger than 0.3 cm in the transverse plane and 0.3 cm 
in the longitudinal direction. For the selected tracks the 
reconstruction efficiency and remaining contamination is 
estimated by Monte Carlo simulations of HIJING [3D] 
and Thcrminator [3U [32] events using a GEANT3 [33] 
detector simulation and event reconstruction. The recon- 
struction efficiency for tracks with 0.2 < p t < 1.0 GeV/c 
increases from 60 to 70% after which it stays constant 
at 70 ± 5%. The contamination from secondary interac- 
tions and photon conversions is less than 5% for p t = 0.2 
GeV/c and less than 1% for p t > 1 GeV/c. Both the 
efficiency and contamination as a function of transverse 
momentum do not change significantly as a function of 
multiplicity and are therefore the same for all centrality 
classes. 

An alternative analysis was performed with tracks re- 
constructed using only the TPC information. For these 
tracks the same selections were applied except for the 
requirement of hits in the ITS and allowing for a larger 
closest distance to the primary vertex, smaller than 3.0 
cm in the transverse plane and 3.0 cm in the longitudinal 
direction. The reconstruction efficiency for these tracks 
with 0.2 < p t < 1.0 GeV/c increases from 70 to 80% after 
which it stays constant at 80 ± 5%. The contamination 
is less than 6% at p t — 0.2 GeV/c and drops below 1% 
at p t > 1 GeV/c. For this track selection the efficiency 
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and contamination as a function of transverse momen- 
tum also do not depend significantly on the track den- 
sity and are therefore the same for all centrality classes. 
The relative momentum resolution for tracks used in this 
analysis was better than 5%, both for the combined ITS- 
TPC and TPC-standalone tracks. The results obtained 
from the ITS-TPC and TPC standalone tracking are in 
excellent agreement. Due to the smaller corrections for 
the azimuthal acceptance, the results obtained using the 
TPC standalone tracks are presented in this Letter. 
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FIG. 2. (color online) a) «2(pt) for the centrality bin 40-50% 
from the 2- and 4-particle cumulant methods for this mea- 
surement and for Au-Au collisions at ys NN = 200 GeV. 
b) U2{4}(pt) for various centralities compared to STAR mea- 
surements. The data points in the 20-30% centrality bin are 
shifted in p t for visibility. 



The ^-differential flow was measured for different 
event centralities using various analysis techniques. In 
this Letter we report results obtained with 2- and 
4-particle cumulant methods denoted v 2 {2} and 

t>2{4}. To calculate multiparticle cumulants we used 
a new fast and exact implementation [33]. The v 2 {2} 
and t>2{4} measurements have different sensitivity to flow 
fluctuations and nonflow effects - which are uncorrelated 
to the initial geometry. Analytical estimates and results 
of simulations show that nonflow contributions to i>2{4} 
are negligible [36) . The contribution from flow fluctua- 
tions is positive for v 2 {2} and negative for v 2 {A} [37] • For 
the integrated elliptic flow we also fit the flow vector dis- 
tribution [38] and use the Lee- Yang Zeroes method [39] . 
which we denote by u 2 {q-dist} and v 2 {LYZ}, respec- 
tively [30] ■ In addition to comparing the 2- and 4-particle 
cumulant results we also estimate the nonflow contribu- 



tion by comparing to correlations of particles of the same 
charge. Charge correlations due to processes contribut- 
ing to nonflow (weak decays, correlations due to jets, etc.) 
lead to stronger correlations between particles of unlike 
charge sign than like charge sign. 

Figure [2^, shows v 2 {p t ) for the centrality class 40-50% 
obtained with different methods. For comparison, we 
present STAR measurements [HI U2] for the same cen- 
trality from Au-Au collisions at ^/s NN = 200 GeV, 
indicated by the shaded area. We find that the value 
of v^ipt) does not change within uncertainties from 
v /s^ = 200 GeV to 2.76 TeV. Figure^ presents v 2 (p t ) 
obtained with the 4-particle cumulant method for three 
different centralities, compared to STAR measurements. 
The transverse momentum dependence is qualitatively 
similar for all three centrality classes. At low p t there is 
agreement of v 2 (p t ) with STAR data within uncertain- 
ties. 
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FIG. 3. (color online) Elliptic flow integrated over the pt 
range 0.2 < p t < 5.0 GeV/c, as a function of event cen- 
trality, for the 2- and 4-particle cumulant methods, a fit of 
the distribution of the flow vector, and the Lee- Yang Zeroes 
method. For the cumulants the measurements are shown for 
all charged particles (full markers) and same charge particles 
(open markers). Data points are shifted for visibility. RHIC 
measurements for Au-Au at ^/s NN = 200 GeV, integrated 
over the p% range 0.15 < pt < 2.0 GeV/c, for the event plane 
U2{EP} and Lee- Yang Zeroes are shown by the solid curves. 



The integrated elliptic flow is calculated for each cen- 
trality class using the measured v 2 (pt) together with the 
charged particle ^-differential yield. For the determi- 
nation of integrated elliptic flow the magnitude of the 
charged particle reconstruction efficiency does not play 
a role. However, the relative change in efficiency as a 
function of transverse momentum does matter. We have 
estimated the correction to the integrated elliptic flow 
based on HIJING and Therminator simulations. Trans- 
verse momentum spectra in HIJING and Therminator 
are different, giving an estimate of the uncertainty in the 
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correction. The correction is about 2% with an uncer- 
tainty of 1%. In addition, the uncertainty due to the 
centrality determination results in a relative uncertainty 
of about 3% on the value of the elliptic flow. 

Figure [3] shows that the integrated elliptic flow in- 
creases from central to peripheral collisions and reaches 
a maximum value in the 50-60% and 40-50% centrality 
class of 0.106 ± 0.001 (stat) ± 0.004 (syst) and 0.087 
± 0.002 (stat) ± 0.003 (syst) for the 2- and 4-particlc 
cumulant method, respectively. It is also seen that the 
measured integrated elliptic flow from the 4-particle cu- 
mulant, from fits of the flow vector distribution, and from 
the Lee- Yang Zeroes method, are in agreement. The 
open markers in Fig. [3] show the results obtained for 
the cumulants using particles of the same charge. The 
4-particle cumulant results agree within uncertainties for 
all charged particles and for the same charge particle data 
sets. The 2-particle cumulant results, as expected due to 
nonflow, depend weakly on the charge combination. The 
difference is most pronounced for the most peripheral and 
central events. 
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FIG. 4. (color online) Integrated elliptic flow at 2.76 TeV 
in Pb-Pb 20-30% centrality class compared with results from 
lower energies taken at similar centralities |40l I43] . 



The integrated elliptic flow measured in the 20-30% 
centrality class is compared to results from lower ener- 
gies in Fig. [4] For the comparison we have corrected the 
integrated elliptic flow for the p t cutoff of 0.2 GeV/c. The 
estimated magnitude of this correction is 12 ± 5% based 
on calculations with Therminator. The figure shows that 
there is a continuous increase in the magnitude of the 
elliptic flow for this centrality region from RHIC to LHC 
energies. In comparison to the elliptic flow measure- 
ments in Au-Au collisions at y/s NN — 200 GeV we 
observe about a 30% increase in the magnitude of vi at 
^/s NN = 2.76 TeV. The increase of about 30% is larger 



than in current ideal hydrodynamic calculations at LHC 
multiplicities [7] but is in agreement with some models 
that include viscous corrections which at the LHC be- 
come less important [121 IT5¥T5] . 

In summary we have presented the first elliptic flow 
measurement at the LHC. The observed similarity at 
RHIC and the LHC of ^-differential elliptic flow at low 
Pt is consistent with predictions of hydrodynamic mod- 
els [3 [2j . We find that the integrated elliptic flow in- 
creases about 30% from ^/s NN = 200 GeV at RHIC to 
y/s NN = 2.76 TeV. The larger integrated elliptic flow 
at the LHC is caused by the increase in the mean p t . 
Future elliptic flow measurements of identified particles 
will clarify the role of radial expansion in the formation 
of elliptic flow. 
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